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Abstract 

In this paper we study generalized classes of volume preserving multidimensional 
integrable systems via Nambu-Poisson mechanics. These integrable systems belong to 
the same class of dispersionless KP type equation. Hence they bear a close resemblance 
to the self dual Einstein equation. Recently Takasaki-Takebe provided the twistor 
construction of dispersionless KP and dToda type equations by using the Gindikin's 
pencil of two forms. In this paper we generalize this twistor construction to our systems. 
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1 Introduction 



In this article we study volume preserving diffeomorphic integrable hierarchy of three flows 
[1]. This is different from the usual two flows cases and this can be studied via Nambu- 
Poisson geometry. It is already known that a group of volume preserving diffeomorphisms 
in three dimension plays a key role in an Einstein-Maxwell theory whose Weyl curvature is 
self-dual and whose Maxwell tensor has an algebraically anti self-dual. Later Takasaki [2] 
explicitly showed how volume preserving diffeomorphisms arises in integrable deformations 
of self- dual gravity. 

Nambu mechanics is a generalization of classical Hamiltonian mechanics, introduced by 
Yoichiro Nambu [3]. At the begining he wanted to formulate a statistical mechanics on 
R 3 , emphasizing that the only feature of Hamiltonian mechanics one should preserve is the 
Liouville theorem. He considered the following equations of motion 

dr 

— = Vti(r)AVt)(r), r = (x,y,z)eR, 

where x, y, z are dynamical variables and u, v are two functions of r. Then Liouville theorem 
follows from the identity 

V • (Vu(r) A Vu(r)) = 0. 
He further observed that the above equation of motion can be cast into 

d£ d(f,u,v) 
dt d(x,y,z)' 

the Jacobian of right hand side can be interpreted as a generalized Poisson bracket. 

Hence the binary operation of Poisson bracket of Hamiltonian mechanics is generalized 
to n-ary operation in Nambu mechanics. Recently Takhtajan [4,5] has formulated its basic 
principles in an invariant geometrical form similar to that of Hamiltonian mechanics. 

In this paper we shall use Nambu mechanics to study generalized volume preserving 
diffeomorphic integrable hierarchy. These classes of integrable systems are closely related to 
the self dual Einstein equation, dispersionless KP equations etc. In fact we obtain a higher 
dimensional analogue of all these systems. It turns out that all these systems can be written 
in the following form: 

d& n) = 

qWaqw = 0. 

For n = 2 we obtain all the self dual Einstein and dKP type equations. 

Hence we obtain a common structure behind all these integrable system, so there is a 
consistent and coherent way to describe all these systems. Here we unify all these classes of 
integrable systems by Gindikin's pencil or bundle of forms and Riemann-Hilbert problem ( 
twistor description). Gindikin introduced these technique to study the geometry of the solu- 
tion of self dual Einstein equations. Later Takasaki- Takebe [6,7] applied it to dispersionless 
KP and Toda equations. 
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This paper is organized as follows: 
In section 2 we give some working definitions of Nambu-Poisson geometry. In section 3 
we construct our volume preserving integrable systems via Nambu-Poisson geometry. If 
one carefully analyse these set of equations then one must admit that they bear a close 
resemblance to the volume preserving KP equation, so far nobody knows about this equation. 
It is known that area preserving KP hierarchy (= dispersionless KP hierarchy) plays an 
important role in topological minimal models ( Landau- Ginzburg description of the A-type 
minimal models ). So we expect volume preserving KP hierarchy may play a big role in low 
dimensional quantum field theories. Section 4 is dedicated to the twistor construction of 
these systems. 

I would like to end this introduction by expressing some regrets. I apologize to the 
mathematicians that the presentation of this paper is not regorous, simply because, I want 
to show the application of the Nambu-Poisson geometry to the large spectrum of readers. 

2 Nambu-Poisson Manifolds 

In this section we will state some basic results of Nambu-Poisson manifold from the paper 
of Takhtajan [4]. 

Let Ai denote a smooth finite dimensional manifold and C°°(M) the algebra of infinitely 
differentiable real valued functions on M.. Recall that [4,5] M. is called a Nambu-Poisson 
manifold if there exits a IR-multi-linear map 

{,..., } : [C°°(M)f n -> C°°(M) (1) 

called a Nambu bracket of order n such that V/i, / 2 , . . . , hn-i £ C°°{M), 

{/!,...,/„} = (-l) £(CT) {/.(l),---,/.(n)}, (2) 
{flh, h, ■ ■ ■ , fn+l} = fl{h, h, ■ ■ ■ , fn+l} + {fl, h, ■ ■ ■ , fn+l}h, (3) 

and 

{{fl, • • • ) fn-1, fn}, fn+l, ■ ■ ■ , hn-l} + {fn, {fl, ■ ■ ■ , fn-1, fn+l}, fn+2, ■ ■ ■ , hn-l} (4) 
H h {fn, ■ ■ ■ , hn-2, {fl, ■ ■ ■ , fn-1, /2n-l}} = {fl, ■ ■ ■ , fn-1, {fn, ■ ■ ■ , hn-l}}, 

where o G S n — the symmetric group of n elements — and e(<x) is its parity. Equations (2) 
and (3) are the standard skew-symmetry and derivation properties found for the ordinary 
(n = 2) Poisson bracket, whereas (4) is a generalization of the Jacobi identity and was called 
in [2] the fundamental identity. When n = 3 this fundamental identity reduces to 



{h, h, h}, U, h} + {h, {h, h, U}, h} + {/ 3 , U, {h, h, h}} (5) 

= {fi, h, {h, h, h}}- 
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It is also shown in [2] that Nambu dynamics on a Nambu-Poisson phase space involves 
n— 1 so-called Nambu-Hamiltonians Hi, ... , H n _i e C°°(.M) and is governed by the following 
equations of motion 

f t ={f,Hi,...,H n _i}, VfeC°°(M). (6) 

A solution to the Nambu-Hamilton equations of motion produces an evolution operator 
Ut which by virtue of the fundamental identity preserves the Nambu bracket structure on 
C°°(M). 

Definition 2.1 / e C°°(Ai) is called an integral of motion for the system if it satisfies 

{f, H 1 ,H 2 , • • • ,H n _i} = 0. 

Like the Poisson bivector the Nambu bracket is geometrically realized by a Nambu 
polyvector 77 e T(A n TM), a section of A n TM, such that 

{fi,---,fn} = v(dfi,---,df n ), (7) 

which in local coordinates (xi, . . . ,x n ) is given by 

, = (x)±- A- A (8) 

where summation over repeated indices are assumed. 

It was stated in [4] that the fundamental identity (4) is equivalent to the following alge- 
braic and differential constraints on the Nambu tensor r\ ix ...i n (x): 

Sij + P(S)ij = 0, (9) 

for all multi-indices % = {ii, . . . , i n } and j = {ji, . . . ,j n } from the set {1, . . . , N}, where 

Sij — Vil-.-inVjl-.-jn VjnW3---inVjl.-.jn-li2 + ' ' ' + Vjni2---in-lhVjl--.jn-lin ~~ Vjni'2-- -inVjl ■■■jn-lh ' (l^) 

and P is the permutation operator which interchanges the indices i\ and j\ of 2n-tensor S, 
and 

V ( n, ■ ' u -'■ +n- . dr lh-jn-ii2 , , .d^i-jn-i^ \ 

iV a 

for all i 2 ,...,i n , ji,...,j n — 1, . . . , TV. 

It was conjectured in [2] that the equation SV,- = is equivalent to the condition that 
n-tensor 77 is decomposable, recently this has been proved by Alekseevsky and author [8], and 
more recently by Marmo et. al. [9] and many others, so that any decomposable element in 
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A n V, where V is an iV-dimensional vector space over R, endows V with the structure of a 
Nambu-Poisson manifold. 

Example 

Let us illustrate how Nambu-Poisson mechanics works in practise. The example is the 
motion of a rigid body with a torque about the major axis introduced by Bloch and Marsden 
[10]. 

Euler's equation for the rigid body with a single torque u about its major axis is given 

by 

thi = aim 2 m 3 (12) 

m 2 = a 2 m 1 m 3 (13) 

m 3 = a 3 mim2 + u (14) 

where u = — km\m<i is the feedback, a\ = — j^, a 2 = — yj- and a 3 = — -jk We 
assume I\ < I 2 < h- 

These equations can be easily recast into generalized Nambu-Hamiltonian equations of 
motion: 

^ = {4^,4 (is) 

where these equations involve two Hamiltonians and these are 

1 a 3 — k 2 2\ 
^2 = t( mi -m 3 ). 



When 

ai = a 2 = a 3 = 1 

and u 3 = 0, these set of equations reduce to a famous Euler equation or Nahm's equation 

(IT- 

= e ij'fe[^})^fe] hjjk= 1,2,3, 

where TjS are SU(2) generators. 

Acknowledgement: Author would like to record here the debt which the work owes 
to the papers of K. Takasaki and T. Takebe and L. Takhtajan and to many illuminating 
discussions with Professors K. Takasaki. He is also grateful to Professors Dmitry Alekseevsky, 
Simon Gindikin, Warner Nahm and Ian Strachan for valuable discussions. 
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3 Volume preserving integrable systems 



In this section we shall follow the approach of Takasaki-Takebe's [6,7] method of studing area 
preserving diffeomorphic ( or sDiff(2) ) K.P. equation. In fact, they adopted their method 
from self dual vacuum Einstein equation theory. In the case of self dual vacuum Einstein 
equation, and hyperKahler geometry also area preserving diffeomorphism appear, where the 
spectral variable is merely parameter. But in case of sDiff(2) K.P. equation the situation 
is different, where Takasaki-Takebe showed that one has to treat A as a true variable and 
it enters into the definition of the Poisson bracket. The situation is similar here, onle the 
Kahler like two form and the associated "Darboux coordinates" is replaced by volume form 
and the Poisson bracket is replaced by its higher order Poisson bracket called Nambu bracket. 

Suppose we consider L = L(\,p, q), M — M(\,p, q) and N = N(\,p, q) are some Laurent 
series in A with coefficients are functions of p and q. 

Definition 3.1 The volume preserving integrable hierarchy is defined by 

dL_ 

8t n 

dM 



dt r< 



= {B ln ,B 2n ,L} (16) 

= {B ln ,B 2n ,M} (17) 
= {B ln ,B 2n ,N} (18) 



dN 
dt n 
and 

{L,M,N} = 1, (19) 

where B in : = (L ra ) n > and B 2n ; = (M n ) n > . The first three equations are hierarchy 
equations and these constitute three flows of the system and the last one shows the volume 
preservation condition. 

Let us now compare our case with the area preserving KP hierarchy. The sdiff(2) KP 
hierarchy which is given by 

^- = {B n ,C} ^ = {B n ,M} {C,M} = 1, 
where £ is a Laurent series in an indeterminant A of the form 

oo 

C = X+J2u n+ l(t)\- n , 
n=l 

B n = (£)>o- The function M. is called Orlov function and it is defined by 

oo oo 

M = ]T nt n C n - 1 +x + Y, ViC^' 1 

71=1 1=1 

where t\ = x. 
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Remark 3.2 Our hierarchy has a structure of volume preserving KP hierarchy and instead 
of one Orlov function we need two Orlove functions Viz. M. and Af. 

Since the flows are commuting, we obtain 

Proposition 3.3 The Lax equation for L, M , N are equivalent to the following equations: 

{ B 2n } + {Bi n , o , 2 " } + {Bi: B 2m } 

ot m ot m 

-{ , B 2m } - {B lm , q^" 1 } + {Bim, H 2 } = 0, 



where 
and 



Hi — {Bin, B 2n , -Blm} 



H 2 — {Bi n , B 2n , B 2m }. 



Proof: Result follows from the compatibility conditions of hierarchy equations and the 
fundamental identity. 
□ 



Remark 3.4 In the case ordinary Poisson geometry and sDiff(2) hierarchy this is a zero 
curvature equation. 

Let Q be a three form given by 

oo oo 

Q :— ^2 dBin A dB 2n A dt n = d\ A dp A dq + dB ln A dB 2n A dt n . 

n=l n=2 

From the definition it is clear is closed 3 form. In fact sDiff(3) structure is clearly ex- 
hibited from this structure and the theory is integrable in the sense of nonlinear graviton 
construction. This is a generalization of nonlinear graviton construction. 

Theorem 3.5 The volume preserving hierarchy is equivalent to the exterior differential equa- 
tion 

Q = dL A dM A dN. 

Proof: We have seen that Q can be written in two ways. Expanding both sides of the 
exterior differential equation as linear combinations of dXAdpAdq, d\AdpAdt n , d\AdqAdt n 
and dp A dq A dt n . 

When we pick up the coefficients of dX A dp A dq, we obtain the volume preserving condition 

{L,M,N} = \. 
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When we equate the other coefficients, viz. d\Adp A dt n , dX A dq A d£ n and dpAdqA dt n 
we obtain the following identities: 

D & \ £>( T 1\/T AT\ 

(20) 
(21) 



and 



d(Bi n , B 2n ) 


d(L,M, N) 


d{\p) 


d(X,p,t n ) 


d(Bi n , B 2n ) 


d(L, M, N) 


d(X,q) 


d(X,q,t n ) 


d(Bi n , B 2n ) 


d(L, M, N) 


d(p, q) 


d(p, q, t n ) 



(22) 



respectively. We multiply the above three equations by ^ and |y respectively. If we 
add first and third equations and substract the second one from them, then after using the 
volume preserving identity we obtain, 

7^3— = {Bin, B 2n , L}. 



Similarly we can obtain the other equations also, in that case we will multiply the equations 

dM dM anr | dM_ 
dq ' dp dUU d\ 



(20), (21) and (22) by f^, and ^ respectively. 



□ 

We can write down the fundamental relation 

= dL A dM A dN 

by 



d(MdL AdN +J2 B Vn dB 2n A dt n ) = 0. (23) 

n=l 

This implies the existence of one form Q such that 

oo 

dQ = Md(LdN) + J2Bi„d(B 2n dt n ). (24) 



n=l 



This is an analogue of "Krichever potential" in the volume preserving case. Hence we 
can say from (24) 

dQ 

M = d(LdN) ' B2 ™>*™ fixed ^ 

Bln = d{B 2 Q n dt n ) ^,N,B 2m ,t m ( m ^n) fixed • ( 26 ) 
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4 Application to mulitidimesional integrable systems 
and Riemann-Hilbert problem 

We already stated that our situation is quite similar to nonlinear graviton construction of 
Penrose [11] for the self dual Einstein equation. This is a generalization of nonlinear graviton 
constructions. 

To the geometer self dual gravity is nothing but Ricci flat Kahler geometry and it is 
characterized by the underlying symmetry groups sDiff(2), this are called area preserv- 
ing diffeomorphism group on surfaces. These are the natural generalization of the groups 
DiffXS 1 ), diffeomorphism of circle. 

It is well known how the area preserving diffeomorphism group appears in the self dual 
gravity equation. Let us give a very rapid description of this. 

Let us start from a complexified metric of the following form 



ds = det 



e 11 e 12 



e 21 e 22 



11 22 12 21 

= e e — e e , 

where e y are independent one forms. Ricci flatness condition boils down to the closedness 
of 

dn kl = (27) 

of the exterior 2-forms 



n kl = -J/Ae 3i , (28) 



1 

2' 

where (J) is the normalized symplectic form, i.e. it is a 2 x 2 matrix whose entries are 0, 1, 
— 1 and respectively. Then above system of two forms can be recast to 



n{\) = ^-J ij (e il + e i2 X)A(e jl + e j2 \). (29) 

This satisfies 

fi(A)AQ(A) = (30) 

dn(\) = 0, (31) 

where d stands for total differentiation. These suggest us to introduce a pair of Darboux 
coordinates 

= dPAdQ 

and these are the sections of the twistor fibration 

7T : T — >CP\ 

where T is the curved twistor space introduced by Penrose. Basically each fibre is endowed 
with a symplectic form and as the base point moves this also deforms and here comes the 
area preservation. 
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Two pairs of Darboux coordiantes are related by 

f(X,P(X),Q(X)) = P' 

g(X,P(X),Q(X)) = Q> 

and / and g satisfy {/, g} = 1. The pair (/, g) is called twistor data. Locally / and g ( 
after twisting with A ) yield patching function. Ricci flat Kahler metric is locally encoded 
in this data. This set up is nothing but the Riemann-Hilbert problem in area preserving 
diffeomorphism case. 

The novelty of this approach is that this twistor construction will work in the higher 
dimensions too, when there is no twistor projection. The most important example is the 
electro-vacuum equation, volume preserving diffeomorphism groups in three dimension play 
a vital role here. This model was first introduced by Flaherty [12] and later Takasaki [2] 
showed how this works explicitly. 



4.1 Gindikin's bundle of forms 

We already stated that anti-self dual vacuum equations govern the behaviour of complex 4- 
metrics of signature (+,+,+, +) whose Ricci curvature is zero and whose Weyl curvature is 
self dual. These two curvatures are independent of change of coordinates, so in one particular 
of the equations these metric becomes autometically Kahler and can be expressed in terms 
of a single scalar function Q, the Kahler potential. Then curvarure conditions will lead you 
to 1st Plebenski's Heavenly equation 

1 (32) 



dxdx dydy dxdy dydx 
and the correspoding anti-self-dual Ricci flat metric is 

d 2 Vt . . . 
g{®) = d iQ j dxW, x' = x,y, x J =x,y. 

This system is completely integrable and it is an example of a multidimensional integrable 
system. 

Let Q be the 2-form, given by 

f2 = dx A dy + X(Vt xx dx A dx + Vt X ydx A dy + VLy^dy A dx + Vt y ydy A dy) + X 2 dx A dy. 



Lemma 4.1 

(1) dtt = 
(2) fiAfi = 0. 
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Proof: Since f2 satisfies Plebenski's Heavenly 1st equation, hence (2) is true. 
□ 



A number of multidimensional integrable systems can be written in terms of a 2-form fl 
which satisfies the equations (4.1). 
Example: 

The dispersionless KP hierarchy has a Lax representation with respect to a series of 
independent ("time") variables t = (ii,i2, • • •) 

f = {B '" C} 

where 

B n :=(C n )> n = 1,2,.-. 
£ is a Laurent series in an indeterminant A of the form 

oo 
n=l 

{ , } is a Poisson bracket in 2D phase space with respect to (A, x). 
Let us consider 

oo 

Q = dX A dx + ^2 dB n A dt n 

n=2 

then fl A Q is equivalent to the zero curvature condition 

dB n dB 



dt m dt 

n 



+ {B n ,B m } = 0. (33) 



This is an alternative form of dispersionless KP hierarchy. 

All these systems are related to area preserving diffeomorphism group sDiff(2). 

In this paper we are presenting an analogous picture for multidimensional integrable 
systems related to volume preserving diffeomorphism group. 

Proposition 4.2 The 3-form 

oo 

fi( 3 ) = dX A dp A dq + dB ln A dB 2n A dt n , 



n=2 



satisfies 



dn (3) = 

Q(3) A Q(3) = 0> 
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This is equivalent to propostion 3.2. Hence it is satisfied. 
□ 

The Gindikin's method [13] of pencil of 2-forms is the most effective way to study these 
systems. Consider the following system of 1st order equation depending on a parameter 
r = (7-1,7-2) GC 2 

e^r) = e 11 r 1 fc H h e lfc r 2 fc , 

e 2 ( r ) = e 21 r k + --- + e 2k r k , 



e 2 \r) = e 2ll r k + --- + e 2lk r k . 
where e* J are 1-forms. Let f2 fc (r) be the bundle of 2-forms 

n k (r) = e 1 (r)Ae 2 (r) + - + e 2,4 (T)Ae 2l (r) 
satisfying the conditions 

.(tt k ) l+1 = 

•(n k y ± 
•dn k = 0. 

The bundle of forms actually encodes the integrability of the original system. In the special 
case I = 1, k = 2, we recover the Ricci flat metric 

g = e ll e 22 _ e 12 e 21 

4.2 Higher dimensional analogue of Gindikin's Pencil 

Let us consider the following system of 1st order equations depending on parameter 

r = (ti, t 2 , 73) G C 2 , 



e 1 ! 


V) 


= e n Tl + e 12 r 2 + e 13 


e 2 ! 


V) 


= e 21 n + e 22 r 2 + e 23 


e 3 i 


V) 


= e 3l n + e 32 r 2 + e 33 



Just like the previous situation these bundle of forms then encodes the integrability of 
the original system. The metric defined here is the higher order analogue of Ricci metric in 
two form case. 

This metric is given by, 



g = e n e 22 e 33 - e n e 32 e 23 + e 12 e 31 e 23 

_ e 12 e 21 e 33 + e 13 e 21 e 32_ e 13 e 31 e 2 2> 
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Remark 4.3 Since 3-form satisfies 



dn® = 

so it denotes the volume preserving multidimensional integrable systems. These systems can 
be described by Gindikin's bundle of multi forms, higher dimensional analogue of nonlinear 
graviton. 

4.3 Twistor description of Volume preserving multidimensional 
integrable systems 

The natural question would be to find out the analogous Riemann-Hilbert problem in the 
volume preserving case. Our situation is very similar to electro-vacuum equation. Let us 
consider two sets of solutions of hierarchy (L, M, N) and (L, M, N) with different analysity. 
Then there exist an invertible functional relation between these two sets of functions such 
that it satisfies 

L = fi(L, M, N), M = f 2 (L,M,N), 
N = f 3 (L,M,N), 

where fi = fi(X,p, q), f 2 = f 2 (X,p, q) and / 3 = /s(A,p, q) are arbitrary holomorphic functions 
defined in a neighbourhood of A = oo except at A = oo. 

We assume fi, f 2 , fz satisfy the canonical Nambu Poisson relation 

{/i,/ 2 ,/ 3 } = I- (34) 

This is a kind of Riemann-Hilbert problem related to three dimensional diffeomorphisms. 
In this case sDiff(3) symmetries is clear, in fact sDiff(3) group acts on (/i, f 2 , ^3), we can lift 
this action on (L, M, N) and (L, M, N) via Riemann-Hilbert fatorization. 

4.4 Application to hydrodynamic type systems 

There are certain kind of integrable systems, called hydrodynamic type, naturally arise in 
gas dynamics, hydrodynamics, chemical kinematics and may other situations, these are given 
by 

u\ = v){u)ul (35) 
where Vj(u) is an arbitrary N x N matrix function of 

u = (m 1 ,---,/), u { = i = 1,---,N. 

The Hamiltonian systems of hydrodynamic type systems considered above have the form 

u\ = K , H}, (36) 
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where H = J h(u)dx, is a functional of hydro dynamic type. The Poisson bracket of these 
systems has the form 

{«*(*) , u»"(y)} = <7 y (ti)«5 x (z - y) + Vi{u)u k J{x - y), (37) 

called Dubrovin-Novikov type Poisson bracket. Dubrovin-Novikov showed if the metric is 
non degenerate i.e. det [<7 y ] 7^ then (37) yields a Poisson bracket provided g^(u) is a 
metric of zero Riemannian Curvature. 

A large class of them can be described by the Lax form, 

ip x = zAip 

ip t = zBip. 

From the compatibility condition we obtain, 

A t = B x 
AB = BA. 
These set of equations can be easily recasted to 

dVL = 
fiAfi = 0, 

where Q is a closed one form. 

Hence has also a twistorial description. The upshot of this section is that Gindikin pencil 
of forms can be applied to large number of classes integrable systems, including the volume 
preserving integrable systems we propose here. 

5 Conclusion 

In this paper we have shown that there are large class of integrable systems can be obtained 
from Nambu-Poisson mechanics. They belong to the same family of self dual Einstein or 
dispersionless KP type equations. In fact they are related to what should be called volume 
preserving KP system. In some sense these integrable systems are higher dimensional gener- 
alization of self dual Einstein equation. Hence these systems are describable via Gindikin's 
bundle of forms, or twistor method. 

In fact we have pointed out in this paper that Nambu-Poisson manifold is an useful tool 
to study volume preserving integrable systems. Recently in membrane theory physicists [14] 
have found M-algebra from M-brane, these are related to Nambu-Poisson mechanics. 

There are certain problems we have not discussed in this paper, viz. the quantization 
of these volume preserving generalized multidimensional integrable systems. Presumably 
method of star product quantization [15,16] would be the best way quantize these systems 
and instead of binary star product we need triple star product [17,18]. Thus it is important 
to see whether Moyal-Nambu bracket can be expressed in terms of Fedosov's triple star 
product. 
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